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ABSTRACT 

We propose a lattice formulation of the chiral fermion which maximally respects 
the gauge symmetry and simultaneously is free of the unwanted species doublers. 
The formulation is based on the lattice fermion propagator and composite oper- 
ators, rather than on the lattice fermion action. The fermionic determinant is 
defined as a functional integral of an expectation value of the gauge current opera- 
tor with respect to the background gauge field: The gauge anomaly is characterized 
as the non-integrability. We perform some perturbative test to confirm the gauge 
covariance and an absence of the doublers. The formulation can be applied rather 
straightforwardly to numerical simulations in the quenched approximation. 
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The chiral fermion on the lattice has refused its manifestly gauge invariant 
treatment [1]. There even exists the no-go theorem [2] for such an endeavor. In 
the continuum counterpart, the chiral fermion develops a curious phenomenon, 
called the quantum anomaly [3] or more definitely the gauge anomaly [4] . It can be 
argued that the difficulty in the lattice chiral gauge theory is a natural consequence 
of the gauge anomaly. 

Suppose that we start with a well regularized fermionic partition function de- 
fined by a manifestly gauge invariant lattice fermion action. In the continuum 
limit, the gauge anomaly is a gauge variation of the partition function. There- 
fore if the fermion content is not free of the gauge anomaly, the partition function 
should not be gauge invariant — this contradicts with the very gauge invariance of 
the formulation. There are two possible resolutions: One is an appearance of the 
species doublers which cancel the gauge anomaly [5]. Another is a pathology in the 
continuum limit such as the non-Lorentz covariance [6]. The "trouble" with the 
lattice regularization is that it always regularizes ultraviolet divergences, even when 
a gauge invariant regularization should be impossible due to the gauge anomaly. 

The above reasoning suggests that the appearance of the unwanted doublers 
is quite natural. However the problem in the conventional approach is of course 
that the doublers appear even in the anomaly free cases. Presumably, the ideal 
lattice formulation of the chiral fermion will be the one which distinguishes the 
anomaly free gauge representations from the anomalous ones. That unknown gauge 
invariant lattice action should have a structure that can be written down, for 
example, for the spinor rep. of so(4) but not for the fundamental rep. of su(3), 
because the latter is anomalous. Such an ideal formulation seems to require a 
further deeper understanding on the origin of the quantum anomaly. 

In this article, we simply abandon a direct gauge invariant definition of the 
fermionic partition function. We take an indirect route. Nevertheless we attempt 
to respect the gauge symmetry as much as possible within a range consistent with 
the gauge anomaly. 
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Instead of directly defining the fermion action and the partition function, we 
start with the propagator and the gauge current operator on the lattice. This 
formulation may be regarded as a first quantization approach, compared to the 
conventional ones. The important fact for us is that although the partition function 
cannot be regularized gauge invariantly in general, the gauge current can always be 
regularized gauge covariantly even if the gauge representation is anomalous. This 
type of regularization scheme in the continuum theory is known as the covariant 
regularization [7]. 

In the covariant regularization, fermion loop diagrams are defined as an ex- 
pectation value of the gauge current J^ a (x), in the presence of the background 
gauge field. The ultraviolet divergence of the diagram is then regularized by in- 
serting a gauge invariant dumping factor into the fermion propagator. In this way, 
the gauge invariance associated with all the gauge vertices except that of J^ a {x) ) 
is preserved. The basic idea is that a possible breaking of the gauge symmetry 
due to the anomaly is forced on the J fxa (x) vertex as much as possible. The gauge 
anomaly (J^ a (x)) thus defined has the covariant form [8,9] because of the gauge 
invariance at external vertices. On the other hand, a gauge singlet operator such 
as the fermion number current is always regularized gauge invariantly. The scheme 
thus spoils the Bose symmetry in general but it is restored when the theory is free 
of the gauge anomaly. The scheme is very powerful and applicable to any chiral 
gauge theories including the Yukawa couplings. 

Once the expectation value of the gauge current is obtained in the covariant 
regularization, the fermionic determinant may be defined as a functional integral of 
(J M °(x)) with respect to the background gauge field. However it is obvious that the 
integration is possible only when there exists a Bose symmetry among all the gauge 
vertices. In other words, the gauge anomaly should satisfy the Wess-Zumino con- 
dition [10] which is a consequence of the integrability. Since the covariant anomaly 
breaks the Bose symmetry and the Wess-Zumino condition [8,9], we cannot define 
the fermionic determinant from the integration of (J^ a (x)), provided that it has 
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the covariant gauge anomaly. Therefore, following the covariant regularization, the 
gauge anomaly is characterized as the non-integrability of this integration process. 
This is a consistent picture because the fermionic determinant cannot be gauge 
invariant when the gauge anomaly is present. Our proposal in this article in spirit 
can be regarded as a lattice version of the covariant regularization. We notice that 
the covariant regularization itself does not require the action level realization [11]. 

We proceed as follows: In the continuum theory, the propagator of a massless 
Dirac fermion is expressed as 

(T1>(x$(y)) = ^r8(x-y) 

lP ! (1) 

where ]J) = 7 M (<9 M +L4 /i ) is the vector type, i.e., non-chiral, covariant derivative and 
F^ v is the field strength. In the second line, the denominator has been rewritten 
as a second derivative, that may allow a lattice propagator free of the doubler's 
massless pole. As the propagator on the lattice, therefore we take 

(T1>(xj$(y)) = G(x,y) 

= iP)(x) - 5(x v) ^ 

~ 1 , n(x) + [^,^][U^x)-l]/(4a 2 ) 

where 5(x,y) = 5 x , y /a and 0(x) is the standard lattice covariant derivative 

m = [U»(x)e ad » - e~ ad »Ul(x)] , (3) 

(Ufj,(x) is the link variable [12] and a is the lattice spacing). To avoid the unwanted 



★ When the gauge group is Abelian, the Wess-Zumino condition is trivial and gives no con- 
straint. The existence of the covariant anomaly implies the non-integrability also in this 
case because, ^ 5d* {J»{x)) /{5A v {y)) ^ d^S {J v {y)) /(SA^x)) = 0, where the left hand 
side is the U(l) gauge anomaly and the right hand side means the current is covariantly 
regularized. 
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massless pole, we define the covariant lattice d'Alembertian by 



U(x) = - -2 \u»(x)e a9 » + e- ad »Ul(x) - 2 



(4) 



For the free theory, this is 2^ (1 — cos ak^)/a 2 in the momentum space and does 



not have the doubler's zero at = n/a. Eq. (4) is nothing but the Wilson term [13] 
apart from one extra 1/a. In eq. (2), U^ v (x) is the standard plaquette variable [12]: 
U^(x) = U^{x)U v {x + a^U^x + a v )ut{x). With the parameterization U^(x) = 
exp[iaAn(x)}, the lattice propagator (2) obviously reduces to the continuum one (1) 
in the naive continuum limit. The choice (2) is by no means unique and other 
definition would work as well. 

Using the lattice propagator (2), we define a fermion bi-linear operator as* 



where the minus sign is due to the Fermi statistics. The gauge current of a right 
handed chiral fermion is simply defined by taking Ai = T a 7 M P# where Pr = 
(1 + 75) /2 is the chirality projection operator. This is possible because nothing 
flips the chirality along the fermion line. When the Yukawa coupling is involved, 
this simple recipe using the Dirac propagator does not work and we will comment 
on later the generalization. 

An important property of the definition (5) is the manifest gauge covari- 
ance. Namely, under the gauge transformation on the link variable U^(x) — > 
V(x)U^(x)V^(x+a^), the propagator (2) is transformed as G(x, y) — > V(x)G(x, y) x 
V*(y). Consequently the bi-linear operator (5) transforms 



which means that the composite operator has a definite transformation property 



* It is equally well easy to define, say, the two point function of baryon type composite 
operators. Note that the present formulation is also applicable to the vector gauge theory 
such as QCD. 




(5) 
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under the gauge transformation on the external gauge field. Note that the covari- 
ance holds for a finite lattice spacing as well as the continuum limit a — > 0. In 
particular, a gauge singlet operator, for which M. commutes with the gauge gener- 
ator, is regularized gauge invariantly. In the continuum limit, the gauge anomaly 
should have the gauge covariant form provided that the limit is not pathological. 

Since we are not assuming the underlying fermion action in the present formu- 
lation, various symmetric properties are unfortunately not manifest. Nevertheless 
the gauge covariance (6) is powerful enough to derive the Ward identity associ- 
ated with external gauge vertices. The vertex function, being a gauge current type 
operator M. = T a ^N inserted, is defined by 



n=l ' j=\ L ^ " • "-■ 



7r la 

4„ 



Xj,fij,aj —n/a 
[■■■^„aa\--a,, 

N 



x 1 \f \PliP2i ■ ,Pn), 



(7) 

where the term independent of the gauge field (n = 0) identically vanishes. For 
example, when the constant matrix M in (7) commutes with the gauge generator, 
we find 

p v Km rjf°»(p) = 0, -ip v lim F^ abc (p, q) - f bad lim I^%) = 0, (8) 

and higher point identities, by examining a variation of the both sides of (7) under 
an infinitesimal gauge transformation. The second relation is nothing but the 
covariant convergence of the gauge current at one of external vertices in the three 
point function. Note however that the gauge covariance (6) itself does not imply, 

<(p M + q,) lrm F% pabc (p, q) - f abd lim T^ dc (p) = 0. (9) 

If this relation would hold, it implies the covariant convergence of arbitrary current 
operators and contradicts with possible anomalies. The crucial point in this formu- 
lation is that vertices associated with the external gauge field and the vertex of the 
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composite current operator are differently treated. Thus in general lim a ^o Tj^ p 
is not symmetric under an exchange /j, «-> v and a <-> b. It is this breaking of 
the Bose symmetry which allows the manifest gauge covariance of the formula- 
tion. However, as was already noted, the Bose symmetry will be restored in the 
continuum limit when the theory is free of the gauge anomaly. 

To verify the above properties and that the unwanted doublers are really absent 
in the present formulation, we explicitly evaluated some of the vertex functions in 
the continuum limit. After a somewhat lengthly calculation using the technique 
in [14], we find that the two point function is given by 



limr^(p) 

log^^- 7 + - + 47T 2 J-— K) , 
—a z p z 3 24 

(10) 

where J = 0.0465 •• • and K = 0.309 •• • are numerical constants [14]. The 
Lorentz covariance is restored and there is no non-local divergence in (10). Also 
the quadratically divergent terms, which are proportional to g^/a 2 , are canceled 
out, as the Ward identity (8) and the hypercubic symmetry indicate. By taking 
Af = Pr, eq. (10) gives the vacuum polarization tensor of a right handed chiral 
fermion: 



'48tt 2 



t v T a T h ^N[jrp v -Yp 2 ) 



1 T Air ^ ^ 

U^ ab (p) = trT a T b (p^ - g^p 2 ) log — ^ - 7 + - + An 2 (J - -K) 

(11) 

It is transverse, as is constrained by the Ward identity (8) and the logarithmic 
divergence has the correct coefficient as a single chiral fermion. Thus we see that 
the formulation in fact respects the gauge covariance and simultaneously is free of 
the species doublers at least in the perturbative treatment. 

For the three point function, we computed the divergence of the vector and 
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the axial gauge currents: 



i(Pn + Qfi) lim T^ pabc (p, q) = 0, 

(12) 

i(p M + 9 M ) Jmnf^g) = ^2 tr T a {T & , T c }e vpal3 p a qp. 

Both relations are consistent with the Ward identity and actually the first of (12) 
may be derived solely from eq. (8) and a general argument. The second relation 
should be interpreted as the covariant anomaly because of the underlying gauge 
covariance: It has the unique covariantized form 

lim (^(x)TVV(z)) = 0, 

Qi >0 

i ( 13 ) 

Therefore the gauge anomaly of a chiral fermion (note that Pr is inserted in the 
gauge current) has the covariant form with the correct coefficient. We also find 
the correct fermion number anomaly [15] of a single chiral fermion by substituting 
T a -> 1 in (13). 

We have observed that, besides the manifest gauge covariance (6), the present 
formulation possesses many desired features at least in the perturbative treatment. 
At this point the reader might be wondering if the present formulation is equivalent 
to a non-local fermion action 



S = a 4 $>(:r) {U(x) + [7*, VHM*) " i]/^ 2 )} (14) 



because it obviously corresponds to the propagator (2), hence the non-locality 
leads to some pathology. This interpretation is not correct. If our formula- 
tion is simply based on the action (14), the gauge current would be defined by 
(ifj(x)T a 'y fI ifj(x))- BS = —5 log j Y\ y dip(y)dijj(y)e s / (5A^(x)) and the definition obvi- 
ously respects the Bose symmetry among all the gauge vertices. As a consequence, 
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we have the consistent form of gauge anomaly in the continuum limit which con- 
tradicts with the manifest gauge invariance of (14). As was already argued, we 
then expect the doublers or, a pathology such as a breaking of the Lorentz covari- 
ance. Our formulation based on the prescription (5) with M. = T a ^^P^ on the 
other hand, explicitly spoils the Bose symmetry but instead respects the maximal 
background gauge covariance. The possible gauge anomaly has the covariant form. 
Therefore two approaches are completely different even in the continuum limit. 

For simplicity of the presentation, we have neglected the possible Yukawa 
couplings up to now, which is important in realistic chiral gauge theories. The 
generalization of (2) and (5) is however straightforward. In the continuum the- 
ory, the covariant derivative is generalized as Tp = 0rPr + PlPl — iG^rPr — 
iGcp^Pi, with obvious notations. The expression (1) is therefore replaced by 
—l/{ip) = ipy (T/ypl) and the following steps are almost identical. The lattice 
d'Alembertian (4) may be used for the right handed and the left handed compo- 
nents respectively. 

From the above analyses, the present proposal seems to provide a gauge co- 
variant (or invariant for a gauge singlet operator) definition of composite operators 
without the unwanted doublers. We then have to integrate the gauge current ex- 
pectation value to construct the fermionic determinant. A cancellation of the gauge 
anomaly is the integrability condition in the continuum limit, as was already noted. 
However this fact is not so useful practically because the analytical integration is a 
formidable task and, the continuum limit is never reached in numerical simulations. 
Clearly we ought to study the integrability with a finite lattice spacing (and asso- 
ciated modifications of (2) and (5), if necessary) for setting up a non-perturbative 
framework. This analysis is in progress and will be reported elsewhere. Here we 
simply note that what is needed in the Metropolis simulation is not the fermionic 
determinant itself but the difference of the determinant between two gauge field 
configurations. This is the lattice analogue of the gauge current expectation value. 

However the integration is not necessary at all if one is contented with the 
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quenched approximation. The application to numerical simulations is straightfor- 
ward once having the lattice fermion propagator such as (2). We therefore believe 
that our proposal, even in the present form, has a range of practical application at 
least within the quenched approximation. 

The author would like to thank K. Fujikawa, S. Kanno and Y. Kikukawa for 
discussions. This work is supported in part by the Ministry of Education Grant- 
in-Aid for Scientific Research, Nos. 08240207, 08640347 and 07304029. 

Note added: After this paper was accepted for publication, the author was 
aware of a similar proposal had already been made [16]. However the point that 
we do not assume the underlying non-local action is the crucial difference. Our 
consideration also explains why the correct axial anomaly evaluated from a com- 
posite current operator definition [16] and the correct vacuum polarization tensor 
are not reproduced in the corresponding non-local action calculations [17]. 
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